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An  Automated  Fingerprint  Identification  System 


J.  H.  Wegstein 

ABSTRACT 

Procedures  are  described  for  automatically  identifying  fingerprints. 
Machine-read  ridge-direction  and  minutiae  data  are  utilized  in  register- 
ing and  enhancing  search  or  file  minutiae  data.     The  quality  of  the 
data  is  measured.     A  procedure  is  then  described  for  utilizing  this 
minutiae  data  in  determining  whether  two  fingerprint  impressions  were 
made  by  the  same  finger. 

Key  Words:     computerized-fingerprint-identification;  identification; 
pattern-recognition . 
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1.  INTRODUCTION 


The  unique  identifying  data  in  inked  fingerprint  impressions  on  standard 
fingerprint  cards  can  now  be  read  automatically  when  the  prints  are  of 
good  quality.-'-'^      The  data  from  latent   (scene  of  crime)  fingerprints 
and  poor  quality  inked  fingerprints  can  be  read  with  semiautomatic 
readers."^      The  first  part  of  this  paper  describes  the  automatic  pro- 
cedure for  preparing  this  data  so  that  it  can  be  stored  in  a  file  or 
searched  against  data  that  has  been  previously  placed  in  the  file.  The 
second  part  of  the  paper  describes  the  automatic  procedure  whereby  a 
computer  can  compare  search  and  file  data  and  determine  whether  or  not 
the  represented  fingerprints  came  from  the  same  person. 

2.     FINGERPRINT  DATA 

The  automatic  fingerprint  readers  at  the  Federal  Bureau  of  Investigation 
can  deliver  both  ridge-direction  and  minute-detail  data  from  each  of 
the  ten  fingerprint  boxes  on  a  fingerprint  card.'^'^  Ridge-direction 
data  consists  of  angles  indicating  the  ridge  directions  at  equally 
spaced  grid  points  superimposed  over  the  entire  print.     The  reader  also 
determines  the  X  and  Y  coordinates  and  the  direction,  9,   for  ridge 
endings  and  bifurcations  (minutiae) .     A  portion  of  a  fingerprint  along 
with  plots  of  ridge  directions  and  minutiae  are  shown  in  Figure  1.* 
Figure  2  shows  how  the  X  and  Y  coordinates,  measured  in  units  of  one- 
tenth  millimeter,  and  the  direction  0,  measured  in  degrees,  are  defined 
for  a  minutia.     If  the  areas  marked  A  are  ridges,  then  the  minutia  is 
a  ridge  ending;     but  if  the  areas  marked  B  are  ridges,  the  pattern  is 
a  bifurcation.     Since  a  ridge  ending  in  one  print  may  appear  as  a 
bifurcation  in  another  print  from  the  same  finger,  no  distinction  is 
made  between  ridge  endings  and  bifurcations  in  recording  data. 

Part  I:     Preparation  of  Fingerprint  Data 
1.  Registration 

Ideally,   the  core  of  a  fingerprint  should  be  in  the  center  of  its 
appropriate  box  and  the  finger  should  be  parallel  to  the  vertical  lines 
of  the  card  when  the  impression  is  made.     In  practice,  the  print  is 
rarely  centered  and  may  be  at  an  angle  of  more  than  20  degrees  with  the 
vertical.     In  order  to  relieve  the  matcher  of  the  burden  of  aligning 
fingerprint  data  at  match  time,  registration  or  centering  and  aligning 
is  performed  on  the  data  before  it  is  placed  in  the  file  or  searched 
against  a  file. 

The  R92  registration  procedure  is  defined  in  Flow  Charts  1  through  6 
and  various  details  will  be  referred  to  by  the  box  numbers  within 
these  charts.     R92  uses  the  29  X  29  array  of  ridge  direction  data, 
ip-L^j,  delivered  by  the  fingerprint  reader  to  find  the  core  or  center 
point   (XX,  YY) .     Figure  3  is  a  plot  of  the  ridge  direction  array  of 
an  ulnar  loop.     The  angles  -90 ._<t|^^90 .  are  in  degrees  and  missing 
angles  have  the  false  value  100  for  easy  recognition  by  the  computer. 

-'-Raised  figures  refer  to  literature  references  appearing  on  page  16. 
^Figures  and  flow  charts  appear  on  pages  17  through  44. 
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The  parameters  used  by  R92  are: 

RK  1  RK  2  RK  3  RK  4  RK  5  RK  6  RK  7  RK  8  RK  9  RK  10  RK  11  RK  12 
111.     45.       0.  8       73.       76.       176.       3       27  1        25  4 


RK  13  RK  14  RK  15  RK  16  RK  17  RK  18  RK  19  RK  20  RK  21  RK  22  RK  23  RK  24 
50.       100.     75.         3  9  3        465.     379.     379.     210.     90.  7. 

RK  25  RK  26  RK  27  RK  28  RK  29  RK  30  RK  31  RK  32  RK  33  RK  34  RK  35 
3  4  3  10        2  76.       45.       30.     180.       0.  2 


1.1     Selection  of  Possible  Center  Points 

Flow  Chart  1  seeks  all  arch-like  patterns  within  the  print  which  might 
house  the  core  or  center  point.     It  begins  at  the  top  row,   i=RK10=l,  and 
sweeps  from  left  to  right  as  j   increases,  seeking  arch-like  patterns 
that  do  not  include  any  blank  spaces.     Boxes  7-15  eliminate  blank 
areas  where  ridge-directions  are  missing.     For  each  area  being  examined: 

i,j-2        i,j-l        i,j  i,j+l      iJ+2  i,j+3 

/       /     /      \      \  \ 

AM  AA  AB  AC  AD  AN 

various  tests  are  performed.     Box  18  requires  that   |  AB  |  + 1  AC[  <_  RKl  (=111°)  . 
Box  20  requires  that  SAD>^RK2   (=45°)  .     Boxes  21-42  accumulate  a 
confidence  score  LCR  which  indicates  if  the  arch  is  well  formed.  For 
example,   if  box  22  finds  that  angle  AA  is  greater  than  0,  box  23 
increases  the  LCR  score  by  1.     If  box  42  finds  that  the  total  LCR  score 
is  less  than  RK4=8,   control  returns  to  boxes  8  and  9  which  moves  one 
step  to  the  right  for  the  next  set  of  ridge  directions. 

If  box  42  finds  LCR  >_  8,   indicating  a  promising  arch,  control  moves 
through  boxes  43  and  44  to  point @  to  Flow  Chart  2.     Boxes  45  to  58 
perform  additional  tests  on  the  ridge  directions  above  and  below 
point   (i,j).     If  a  crossover   {       V^)   is  found  above  in  row  i-1  or  beJov 
in  row  i+1  near  column  j,   then  SW>0,  and  control  goes  to  box  59 
where  the  reference  point  (i,j)  and  other  information  are  entered  in 
the  K  tables.     For  example: 

*TI^  =  i,     TJ^  =  j ,  and 

TSAD  =SAD=|v.    .   J  +  k.    -l  +  k.    .^il  +  k.  .^ol 
K  i,J-l         1,3         1,3+1  1,3+2'. 

In  Figure  4,   showing  a  portion  of  the  ridge  direction  data,  SAD=246. 
Control  then  passes  to  box  59.01  and  on  to  box  78.     During  this  phase 
the  sum  of  the  absolute  values  of  four  consecutive  ridge  directions 
from  row  i-f-l  is  obtained  as  sho\sm  in  Figure  4. 
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STX=|*i+,,j.3l  +  ln+i,j.2l  +  l*i«,j-il  +  ll'tW,j 

Then  j  is  increased  by  1  and  another  STX  is  computed.     This  is  repeated 
until  there  are  five  different  values  for  STX.     The  largest  of  these 
values  is  called  SXH  and  the  lowest  is  called  SXL.     In  Figure  4,  the 
largest  STX  =  335.     Hence  SXH  =  335.     The  lowest  STX  gives  SXL  =  307. 
SDY  =  SXH  -  SAD  =  89.     These  will  be  entered  in  the  K  tables  as 
TSXH^,  TSXL^,  and  TSDY^. 

With  the  above  data,  box  74  finds  that  SXL  =  307  is  greater  than 
RK6  =  76  and  control  goes  to  box  75.     Here  SDY  =  89  is  found  to  be 
larger  than  RK5  ="  73  and  box  76  finds  SXH  =  335  to  be  larger  than 
RK7  =  176.     Thereupon  box  77  sets  TRY^  =  1  indicating  that  this  area 
may  embrace  the  core  of  the  fingerprint  and  that  the  pattern  may  be  a 
loop  or  a  whorl. 

Control  then  returns  via  point  (^to  Flow  Chart  1  where  i  and  j  are 
advanced  so  that  the  entire  print  is  explored  for  all  possible  core 
points  which  are  then  added  to  the  K  table.     After  the  lowest  row  of 
ridge  directions,  as  designated  by  ILOW  and  lYLB,  has  been  tested,  KM 
is  set  to  K  the  total  number  of  entries  in  the  K  table  and  Control 
goes  via  point  (a) to  Flow  Chart  3.     At  this  time  the  K  table  correspond" 
ing  to  the  ridge  directions  shown  in  Figures  3  and  4  is  shown  in  Figure 
5. 

1.2    Analysis  of  Area  Around  Center  Point 

In  Flow  Chart  3,  boxes  84,  85,  and  86  seek  the  last  entry  in  the  K 
table  where  TRY  =1.     In  the  above  example  this  occurs  when  K=8  at  the 
reference  point   (lYP , JXP) = (12 , 16) .     Box  89  then  defines  two  areas 
around  this  reference  point  shown  as  areas  A  and  B  in  Figures  3  and  6. 
Further  tests  are  performed  on  the  ridge-directions  in  these  areas 
working  upward  from  the  bottom  of  box  B.     Boxes  94-131  and  144-145  lead 
to  RAS  in  box  102  which  is  ten  times  the  mean  of  the  absolute  values  of 
the  second  differences  of  the  angles  in  area  B  of  Figure  6. 

Boxes  132-143  deal  with  the  ridge  directions  in  area  A.     IC  in  box  133 
is  the  total  number  of  cases  where  the  absolute  values  of  the  angles 
are  greater  than  RK30(=76?) .   ITG  in  box  141  is  the  total  number  of 
times  that  angles  are  too  large  for  an  arch  form  in  area  A.     IIF  in 
box  143  counts  the  number  of  times  that  the  absolute  values  of  the 
angles  in  any  row  of  area  A  are  less  than  RK24(7°).     In  other  words, 
this  indicates  excessive  flatness  in  any  row.     Boxes  144-154  are  again 
used  to  accumulate  second  differences,  this  time  in  area  A.  Boxes 
105.1  and  105.2  determine  IF,  the  largest  of  the  IIF  values  from  the 
three  rows  in  area  A.     When  box  105.21  finds  that  the  analysis  of  the 
top  row  in  area  A  has  been  completed,  control  goes  via  point^to 
Flow  Chart  4. 
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In  box  110,  RAT  is  ten  times  the  mean  of  absolute  second  differences 
for  area  A.     This  is  similar  to  RAS  for  area  B  in  Figure  6.     The  various 
quantities  computed  in  Flow  Chart  3  are  renamed  and  stored  in  Table  K 
for  possible  later  use.     Boxes  114-124  perform  various  tests  on  these 
quantities  which  may  lead  to  setting  SW3=1  in  box  125  indicating  that 
this  entry  in  the  K  table  is  unsatisfactory  as  a  site  for  the  core.  For 
example,  RASK>RK21,  where  RK21=379,  indicates  a  roughness  in  the  ridge 
directions  like  the  area  around  (4,24)   in  Figure  3  which  would  be 
unacceptable  for  a  core  area.     TIF>_RK26   (4)  would  indicate  that  there 
are  too  many  flat  ridge  directions  inarea  A.     If  any  of  these  tests 
set  SW3=1,  control  upon  arriving  at  box  127  will  pass  via  point   @  to 
box  86  in  Flow  Chart  3.     Then  control  will  back  up  in  the  K  table  and 
try  the  next  entry  where  TRY  =1. 

1.3    Center  for  a  Loop,  Whorl,  or  Tent  (CLU=1) 

In  the  example  being  followed  with  Figure  6,  control  arrives  at  box  127 
with  SW3=0  and  passes  through  point©  to  box  197  in  Flow  Chart  5.  Since 
K=KM,  box  205  sets  CLU=1  indicating  the  pattern  is  probably  a  loop, 
whorl,  or  tent  and  control  passes  via to  Flow  Chart  6  where  the 
precise  center  point   (XX, YY)   is  computed  as  described  in  Section  1.6, 
In  box  197,  K+2j<KM  would  indicate  that  there  is  some  arch  pattern 
below  the  chosen  core  area  and  boxes  198-203  would  set  out  to  verify 
that  this  is  merely  an  arch  pattern  under  a  delta  or  other  harmless 
anomaly  and  return  to  box  205.     Failing  this,  control  would  pass  via 
point       to  box  155  in  Flow  Chart  4  for  examination  as  an  arch  pattern. 

1.4     Center  for  an  Arch  (CLU=2) 

Before  concluding  with  the  computation  of  (XX, YY)  let  us  return  to 
box  85  in  Flow  Chart  3.     If  there  were  no  TRY  =1  in  the  K  table, 
control  would  loop  through  boxes  86  and  84  till  K=0  and  then  pass  via 
point       to  box  155  in  Flow  Chart  4  to  locate  a  center  point  for  an 
arch  pattern.     Figure     7  shows  the  K  table  for  the  portion  of  arch 
pattern  shown  in  Figure  8.     If  box  155  finds  at  least  two  entries  in 
the  K  table,  boxes  157-163  will  start  with  the  reference  point 
(TI2,  TJ2)=(12,18)  and  determine  that  either  the  previous  or  following 
entry  in  the  K  table  will  contain  another  reference  point  that  is  less 
than  two  spaces  away  in  the  ridge  direction  matrix.     This  is  true  for 
all  of  the  (TI,  TJ)   in  Figure  8.     Consequently,  for  all  entries  except 
K=7  and  8  which  are  rejected  by  box  162  for  having  TSXL    less  than  50  , 
control  moves  to  boxes  167-170.     Here  ^  1 |   for  each  row  in  Figure  8  is 
computed  and  stored  as  TSA    in  table  K  shown  in  Figure  7.     For  each 
of  these  entries  that  passH^oxes  172.1  and  172.2  to  point  control 
returns  to  Flow  Chart  3  where  the  neighboring  ridge  directions  are 
examined  as  described  in  Section   1.2.    Note  that  now  SW1=1  so  that 
when  control  successfully  reaches  box  126  in  Flow  Chart  4,  control 
passes  through  point (d) and  on  to  boxes  174  and  175.     Here  TRF  is  set 
to  the  largest  TSA^(=392)   in  the  K  table  and  KT=5  notes  the  entry  in 
the  K  table  of  the  promising  center  point  for  the  arch  pattern  in 
Figure  8.     When  the  last  entry  in  the  K  table  has  b^en  examined. 


box  176  sends  control  via  point ^  to  Flow  Chart  5  for  a  further 
examination  of  the  ridge  directions  under  the  proposed  center  point 
(15,  18)   in  Figure  8.     Boxes  180-194  define  an  area  bounded  by  the 
points  (16,  16),  (16,  21),    (24,  16),    (24,21)  and  count  the  number  of 
times,  JCwhen   \^ .    .|^75o.   If  JC£3 ,  box  195  sets  CLU=2  and  names  the 
center  reference  point   (lYP,  JXP)  =  (TI     ,  TJ^^^)  .     In  the  example  of 
Figure  8,    (lYP,  JXP)  =  (15,  18).     Control  then  goes  via  point  ©  to  Flow 
Chart  6  where  the  precise  center  point   (XX, YY)  is  computed  as  described 
in  Section  1.6» 

1.5     No  Definite  Center  Found  (CLU=3) 

If  boxes  206  and  207  in  Flow  Chart  5  find  CLU>2  and  CLX>0  indicating 
that  there  was  at  least  one  TRY^=1  in  the  K  table  even  though  none  were 

accepted,  then  box  209  sets  KT=KTX  and   (JXP,   lYP)  becomes  a  "guess" 
at  the  center  point.     (In  box  116  in  Flow  Chart  4  KTX  had  noted  the 
last  entry  in  the  K  table  where  TRYj^=l.)     If  CLX=0  indicating  that  there 
was  no  TRY^=1  in  the  K  table,  then  boxes  208.1-208.5  make  one  last 
desperate  effort  to  name  a  center  point  by  selecting  the  last  entry  in 
the  K  table  where  TSDY>36°.     If  this  fails,   (XX,  YY)   is  left  at  (200, 
200)  and  control  passes  from  point  0  to  point  (g),  the  end  of  the  program. 
A  few  numbers  that  might  be  useful  for  diagnostic  purposes  are 
identified. 

1.6     Computation  of  the  Center  Point  (XX,  YY) 

When  control  arrives  at  point  (h) in  Flow  Chart  6,  CLU,  KT,   lYP,  JXP 
been  assigned  values.     The  next  task  is  to  obtain  a  more  precise 
location  for  the  center  or  core  point  in  rectangular  coordinates 
(XX,  YY) .     This  center  point  marked  +  for  three  ridge  patterns  is 
shown  in  Figure  9.     The  rectangular  coordinate  system  that  is  used 
dealing  with  minutiae  is  related  to  the  grid  point  system  of  ridge 
directions  as  follows: 

X(j)  =  .4849(49+24(j-l)) 
Y(j)  =  .4849(725-24(1-1)) 

Boxes  211  and  212  first  find  the  abscissa, 

XXI  =  X+X(j), 

of  a  point  between  ii     .  and  x.    .  ,i   as  shown  in  Figure  10.     For  non- 

i»J  i»J+l 
arches,  noted  by  box  213,  box  219  uses  the  ratio' 

PHI  =  DSPl 
DS  90 

to  obtain  a  measure  of  the  verticalness ,  DHl ,  of  the  ridge  directions 
at   (i,  j)  and   (i,  j+1) .     Similarly,  boxes  223-226  computer  XX2  and  DH2 
for  the  next  lower  row  of  ridge  directions.    (The  crossover  reference 
point  JL=TJL      was  obtained  in  box  59.01.)     Using  DH,  the  mean  of  DHl 
and  DH2 ,  box  227  uses  the  ratio  shown  in  Figure  11  to  obtain  XX. 
Box  228  follows,  obtaining  YY=Y (i) -DS+DH . 
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For  the  loop  pattern  shown  in  Figures  3,  4,  and  6  where  (i,  j)=  (12,16): 

DSP1=43+67=110 
XXI  =  11-6376x43  ^  193.3=202.8 

Since  DSPl >  90,  DH1=DS=11 . 6376 .     Since  there  are  no  negative  ridge 
directions  in  row  13,  JL=0,  XX2=XX1,  DH2=11.6376,  DH=11.6376,  and 
XX=XX2=XX1=202 . 8 .     YY=223 . 5-11 . 6376+11 . 6376=223 . 5 . 


With  an  arch  pattern , where  CLU=2,  control  normally  goes  from  box  212  to 
boxes  215,   215.1,  and  216.1.     Using  the  data  in  Figure  7  which  is 
associated  with  the  arch  in  Figure  8: 


(i,j)=(15,18)  DSP1=16+15=31. 
XX=XX1  =  ^^-^^^^^^^  +  221.6=227.6 


YY  is  obtained  by  a  "center  of  gravity"  approach  as  illustrated  in 
Figure  12: 


TNUM=177x346+189x392+200x366=208530. 


TDNM=346+392+366=1104 


The  registration  concludes  in  Flow  Chart  5  by  assembling  some  additional 
data  for  possible  future  use. 

2.     (MANGA)  Procedure  for  Improving  Minutiae  Angles 

Due  to  flaws  in  the  fingerprint  image,  minutiae  direction  angles  are 
sometimes  in  error  with  respect  to  neighboring  ridge  direction 
angles.     For  example,  see  Figure  13.     Flow  Charts  7  and  8  describe  a 
procedure  for  correcting  these  angles  by  utilizing  the  angles  of  the 
neighboring  ridge  directions.     Box  2  in  Flow  Chart  7  defines  a 
2xDSP  by  2xDSP  square  centered  about  each  minutiae.     Boxes  3-25 
construct  a  table  of  the  up  to  K=4  ridge  directions  (XK^, YK^,AK^)  = 

^'^jl2'  ^j22'  '^j32'^  that  surround  each  minutiae  as  shown  in  Figure  14. 

(Note  that   (A.-,-.,  A..,,  A._,)  =  (X.,  Y  ,  0.)  represents  the  coordinates 

111        izl        iJl  111 

and  angle  of  minutia  i  and   (^^-^2^  ^^^22'  '''^j32'^  represents  the 
coordinates  and  angle  of  ridge  direction  j.)     The  A    _  angles  are 

between  -90°  and  +90°.     The  search  for  ridge  directions  surrounding  a 
minutia  seesaws  back  and  forth  across  the  ridge  direction  table  from 
j=l  to  j=IMX2. 
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Control  now  passes  via  point (b) to  Flow  Chart  8.     In  box  28  the  distance 
between  the  minutiae  and  each  ridge  direction,  as  shown  by  broken  lines 
in  Figure  14,   is  subtracted  from  DAG,  the  diagonal  distance  of  the  DS 
square,  and  set  equal  to  ATj^.     SN  is  the  sum  of  these  weighting  factors 
multiplied  by  the  absolute  value  of  their  corresponding  angles  and  AMN 
is  the  mean  of  these  products.     If  AMN  is  less  than  45°  as  determined 
by  box  33,  the  neighboring  ridge  directions  are  predominately  horizontal 
as  in  Figure  15,  and  boxes  34-39  compute  the  weighted  mean  angle  TP  in 
a  straightforward  manner.     If  AMN  is  greater  than  45°  indicating 
predominately  vertical  ridge  directions  as  in  Figure  16,  negative  angles 
indicated  by  box  40  must  be  converted  to  positive  values  before  arriving 
at  the  weighted  mean  in  box  44.     If  the  absolute  value  of  DA,  the 
difference  between  TP  and  the  minutia  angle,  is  greater  than  DAK(=60°) , 
this  minutia  is  deleted  from  the  list  of  minutiae  and  control  returns 
via  (a)  to  examine  the  next  minutia.     If   |da|       DAK  the  angle  value  TP 
replaces  the  angle  in  the  minutia  table  and  control  returns  via^a) 
to  consider  the  next  minutia. 

3.     (NU14)  Computation  of  Rotation  Angle 

After  the  core  or  center  point   (XX,  YY)  of  a  fingerprint  has  been 

determined,   the  minutiae  data  will  be  rotated  about  this  point  through 

some  angle  V  to  a  standard  orientation  for  filing  or  searching.  The 

computation  of  A/  is  described  in  Flow  Chart  9  and  uses  the  following 
parameters  : 

NUMX    KJC     KJL     KUP  KDN     KCT  KNUl 

20.       2        9        3  6        0  20. 

Box  2  determines  two  areas  about  the  reference  point   (lYP,  JXP)  shown 
as  area  C  and  D  in  Figure  3  and  as  areas  A  and  B  in  Figure  17  using 
the  parameters  KJL,  KUP,  and  KDN.     Within  this  area  boxes  9-13 
systematically  locate  pairs  of  ridge  directions  iIi .    .  and  iIi .    ..  on  the 

same  row  and  equidistant  from  the  center  line  as  seen  in  Figure  17. 
Box  14  computes  RN  which  is  half  the  difference  of  each  of  these  pairs 
of  angles,  and  box  15  saves  only  those  RN  where   |rn|    <30°.     Wnen  all  of 
these  pairs  of  angles  in  areas  A  and  B  of  Figure  17  have  been  used,  'V 
the  mean  value  of  these  RN,   is  computed  in  box  22  followed  by  the  sine 
and  cosine  in  box  24. 

4.     (TRARO)  Translation  and  Rotation  of  Minutiae  Data 

Flow  Chart  10  describes  the  procedure  for  translating  the  minutiae  data 
centered  at   (XX,  YY)   to  the  standard  center  location  (200,  200)  and 
at  the  same  time  rotating  the  minutiae  pattern  through  the  angle  V. 

5.     (CLIP  3)  Elimination  of  Excess  Minutiae 

Figure  19  shows  the  162  raw  minutiae  as  they  came  from  the  fingerprint 
reader.     These  correspond  to  the  ridge  directions  shown  in  Figure  3. 
Procedure  CLIP  3  retains  only  those  minutiae  within  the  area  described 
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in  box  1  of  Flow  Chart  11.     The  results  of  this  clipping  are  shown  in 
Figures  19  and  20}  normally  PXA=80,  PYU=70,  PYD=110. 

6.     (YSORT)  Ordering  Minutiae  Data 

The  high  speed  fingerprint  matchers  M40  and  M82  require  that  both 
search  and  file  minutiae  be  arranged  in  descending  order  of  the  Y  values. 
This  is  accomplished  before  searching  or  filing  the  data  by  the  proce- 
dure described  in  Flow  Chart  12. 


7.     Function  DA(TH1,  TH2) 


The  FORTRAN  expression  DA(TH1,  TH2)  calls  for  the  procedure  described 
in  Flow  Chart  13  which  computes  the  difference,  DA,  between  the  ridge 
direction  angles  THl  and  TH2  making  sure  that  |da|  <_  90°.  It  is  used 
by  the  Q6  procedure. 

8.     (Q6)  Measure  of  Fingerprint  Data  Quality 

The  procedure  described  in  Flow  Charts  14,  15,  and  16  examines  the 
ridge  direction  data  for  poor  quality  in  the  fingerprint.     In  Flow 
Chart  15  boxes  1-19  define  a  154  point  area  in  the  ridge  direction  array 
centered  about  the  reference  point   (lYP,  JXP) .     Then  for  each  point 
(i,  j)   in  this  area,  subroutine  QUAL  shown  in  Flow  Chart  14  is  called 
upon.     In  box  7  the  Qn  differences  between  the  angle  at   (i,  j)  and  the 
angle  above,  to  the  left  and  right  and  below  are  computed  by  function 
DA  described  in  Flow  Chart  13.     (Ridge  direction  'B^^   is  the  same  as 

i|j .  .  used  in  earlier  flow  charts.)     These  Qn  are  first  differences 

comparable  to  first  partial  derivatives  — ^  and-rrr.     Box  8  obtains  QA, 

OA  0  1 

the  mean  of  the  absolute  values  of  these  Qn.     In  box  9  the  second 

3  B  3  ^B 

differences  QX  and  QY,  comparable  to    -^^7  ^'^^  'y^I'  obtained.  The 

mean  of  the  absolutes  of  these,  Q2A,  along  with  QA  are  the  outputs  of 
this  subroutine.     MQ=1  indicates  that  both  values  could  be  computed. 
The  results  for  two  ridge  configurations  about   (i,  j)  are  shown  in 
Figure  18.     A  large  value  for  Q2A  indicates  a  disturbed  ridge  direction 
area.     Boxes  15  and  21  of  Flow  Chart  15  obtain  the  mean  values  QAM  and 
Q2AM  for  the  entire  area  of  ridge  directions  defined  earlier. 

The  second  part  of  Q6  in  Flow  Chart  16  uses  the  following  parameters: 

KQl      QX2      QK3      QK4      QK5      KQ6      KQ7      KQ8      KQ9      KQIO  KQll 
120      17.       19.       20.5     72.        5  95        80        40        20  15 

Working  with  the  ridge  directions  of  each  fingerprint,  it  produces  a 
number  of  demerits,  ND,  ranging  from  0  to  3 .     If  either  of  a  pair  of 
mating  fingerprints  that  are  being  matched  by  computer  has  an  ND >  0 , 
the  resulting  matching  score  will  tend  to  be  low.     The  larger  ND 
becomes,   the  poorer  will  be  the  matching  score.     An  ND=3  predicts 
failure  to  identify  the  fingerprint  and  suggests  that  the  minutia 
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should  be  reread  with  a  semlautomated  reader. 

Boxes  24-26.4  accumulate  a  number  of  demerits  depending  on  how  many  of 
the  ridge  directions  in  the  154  areas  being  examined  by  QUAL  are  missing. 
Boxes  27-32  begin  accumulating  demerits  when  Q2AM  is  greater  than  16. 
If  both  XX  and  YY  are  200.,  ND  is  immediately  set  to  3  by  box  35.  Boxes 
41-49  count  the  number  of  cases,  IV,  where  the  absolute  values  of  the 
ridge  directions  exceed  72°  in  the  narrow  area  above  the  center  point 
defined  by  box  38.     If  box  50  finds  IV  to  be  five  or  more,  box  51  sets 
ND=3.     Finally,  boxes  52-57  give  demerits  if  the  total  number  of 
minutiae,   IMG,   is  low.     In  conclusion,  boxes  58  and  59  limit  the  total 
number  of  demerits  for  a  fingerprint  to  three. 

9.     (PRO)  Profusion  of  Minutia 

Subroutine  PRO  described  in  Flow  Chart  17  constructs  a  square  box 
centered  at  each  minutia.     The  box  measures  40x40  units  when  parameter 
DK=20.     The  procedure  then  counts  all  of  the  MT  minutia  in  this  box. 
After  all  of  the  minutiae  in  the  fingerprint  have  been  thus  examined, 
MD  is  set  equal  to  the  largest  MT  that  was  found.     Thus,  MD  is  a 
measure  of  the  maximum  minutiae  density  found  in  the  print. 

'   .■  10.     (TIUP)  Minutiae  Tidy  Up 

When  two  minutiae  are  too  close  together,  one  or  both  minutiae  are 
eliminated.       First,  boxes  1-3  in  Flow  Chart  18  clear  an  array  L.  Its 
length  is  IMO  which  is  equal  to  the  number  of  minutiae  being  examined. 
Boxes  5-22  describe  the  procedure  for  detecting  two  minutiae,  i  and  K, 
which  are  less  than  TK(=10)units  apart  in  X  and  Y.     If  this  occurs  and 
both  minutiae  are  pointing  within  TA(~30*-')  of  each  other,  a  note  is 
made  in  array  L  by  setting  Lj^=l  in  box  17.     If  the  pair  of  minutiae  are 

pointing  in  opposite  directions,  both  L.  and  L    are  set  equal  to  one  in 

1  K 

box  20.     Boxes  26-31  then  eliminate  all  minutiae  identified  by  a  one  in 
array  L. 

11.     (TIUP  1)  Another  Version  of  TIUP 

TIUP  1  is  the  same  as  TIUP  except  that  box  20  contains  only  L     , .  In 

K=i 

the  case  of  opposing  minutiae,  only  the  second  one  is  eliminated  and 
this  appears  to  produce  slightly  better  matching  scores  for  mating 
prints . 

A  comparison  of  Figures  19  and  20  shows  the  results  of  all  of  the 
preparation  procedures  concluding  with  TIUP  1.     For  example,  minutiae 
104  and  105  were  too  close  together  in  Figure  19.  Consequently, 
minutia  105  was  deleted  and  only  104  remains  as  minutia  51  in  Figure  20. 
The  X,  Y,  e  data  corresponding  to  Figure  20  can  now  be  placed  in  the 
fingerprint  file. 

This  concludes  the  preparation  of  minutiae  for  searching  or  filing 
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unless  yet-to-be  defined  procedures  are  developed  for  classification 
and  compaction  of  the  minutiae  data. 


Part  II.     The  M82  Fingerprint  Matcher 

1.     The  MAO  Matcher 

The  M82  fingerprint  matcher  is  an  extension  of  the  MAO  matcher  reported 
in  an  earlier  publication^.     This  extension  utilizes  a  tensor  that 
produces  improved  matcher  scores  particularly  in  cases  where  the  mating 
prints  are  distorted.     Parts  1  and  2  of  the  M82  matcher  shown  in  Flow 
Charts  20  and  21  are  essentially  the  same  as  the  MAO  matcher. 

Figure  21  shows  four  search  minutiae  superimposed  on  five  file  minutiae 
The  file  set  is  displaced,  rotated,  and  stretched  with  respect  to  the 
search  set.     The  X,  Y,  0  data  for  these  minutiae  are  shown  in  Figure  22 
As  M82  sets  out  to  compare  these  prints,  boxes  1-3  in  Flow  Chart  20 
examine  the  total  number  of  search  and  file  minutiae  and  assign  the 
subscript  B  to  the  set  with  the  least  minutiae  (search  set)  and  C  to 
the  other  set   (file  set).     Boxes  5-18.1  proceed  to  put  a  box  around 
each  search  minutia  i  as  shown  in  Figure  23.     For  each  file  minutia  j 
in  this  box  with  V9<LQ  an  entry  S  is  made  in  the  S  table  shown  in 
Figure  2A.     The     indices  i,  j  for  each  pair  of  prints  is  also  entered 
as  ISg  AND  JSg  in  this  S  Table.     (These  values  were  not  saved  in  the 
MAO  matcher.)     The  parameters  used  in  this  illustration  are  as  follows: 

LG      LS       KG      AKR      KHG      KHR      KQL      KMX      XYK      ANK  BKR 
15      25      0        26        16        26        -9  7  5        15.  5 

KRR      PS       BB       CC      DD       EE  KMM 
8500    7        20       .2      30      .36  15 

2.     Computation  of  R  and  TS  Scores 

After  the  last  search  minutiae  i  has  been  examined,  control  passes  via 
(B)to  Flow  Chart  21  where  a  score  R  is  computed.     Boxes  23-23.5  begin 
by  setting  all  values  of  table  TS  to  zero.     Each  entry  i  in  Figure  2A 
is  then  compared  with  each  other  entry  j  and  the  differences  TR  and  T0 
defined  in  box  25  are  computed.     Those  values  of  TR  which  pass  the  test 
of  boxes  26-28  are  used  to  compute  a  score  KR-TR  in  box  28.1    This  is 
added  to  the  TS ,  and  TS .  at  the  locations  corresponding  to  the  minutiae 
pairs  IS  and  JS'^that  contributed  to  this  score.     For  example,  using 
entries  i=l  and  j=2  as  shown  in  box  25: 

TR= 1 1+7  I +  [6-8  I =8+2  =  10 
KR-TR  =  26-10  =  16 
(Note  that  in  box  A  in  Flow  Chart  20,  KR=AKR=26.) 

Then  in  box  26:  TS^=0+16=16 


TS2=0+16=16 
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Stated  in  another  way,   search  minutia  1  in  figure  21  sitting  on  file 
minutia  1  coupled  with  search  minutia  2  sitting  on  file  minutia  2 
produces  a  score  of  16.     This  is  credited  to  the  TS  score  for  each  pair 
in  Figure  24. 


The  same  accumulated  score  called  R  in  box  29  is  the  R  score  used  to 
produce  the  final  RS  score  in  the  M40  matcher. 


3.     Ordering  and  Clipping  the  S  Table 


When  all  combinations  of  the  entries  in  the  S  table  have  been  processed 
control  moves  via  point  (g) to  Flow  Chart  22.     Since  KS¥  was  set  equal  to 
zero  in  box  4  of  Flow  Chart  20,  control  passes  to  box  34.02  where  the 
switch  KSW  is  set  to  1.     With  mating  cards  the  S     table  will  usually 
contain  more  than  100  entries.     It  is  the  task  of  boxes  40-57.1  to 
arrange  the  TS  values  of  the  S  table  in  descending  order  and  retain  the 
associated  IS  and  JS  values.     Thus  IS,  JS,  TS  become  GI,  GJ,  GT  in  the 
J  table  shown  in  Figure  25  with  the  GT  now  in  descending  order.  This 
is  accomplished  by  comparing  each  TS  value  with  each  GT.     When  TS  is 
greater  than  GT ,   this  GT  and  the  remaining  entries  in  the  GT  table  are 
pushed  down  one  notch.     The  TS  value  is  then  inserted  in  the  vacated 
GT  space  and  the  process  is  repeated  till  all  TS  values  are  treated. 
In  practice  only  KMM=15  entries  are  retained  in  the  J  table.     These  are 
normally  mating  minutiae  and  their  TS  contributions  to  the  R  score. 


■  4.     Computation  of  the  Tensor 

Control  next  moves  via @ to  Flow  Chart  23  where  table  J  is  utilized  in 
computing  the  four  element  tensor:     GOV,  SIV,  SIW,  COW. 

The  following  parameters  are  used  in  practical  applications: 

L9  IS         KG  AKR        KH9  KHR      KQL     KMX  XYK        ANK  BKR 

15.47     25.00      .00      8.00      16.00      8.00      -9        7     5.00    15.47  5. 


KRR  PS       AA     BB        CC        DD        EE        BCMM  FRC 

8500.00      7.       0.     20.     .130      40.     .360      15  .0 


For  the  following  illustration,   some  of  these  will  be  changed  as 
follows: 

AKR      KHR      CC  DD 
26.       26.       .20  30. 

Box  58.002  assigns  the  coordinate  labels  (XXU,  YYU)  and  (XXF,  YYF)  to 
the  search  and  file  m.inutiae  that  produced  the  highest  TS  score  by 
using  the  first  entry  in  the  J  table  in  Figure  25.     Using  the  first  GI 
and  GJ  values  and  consulting  Figure  22  one  sees  that  the  search 
minutia  number  2  coordinates  are  (XXU,  YYU)= (181 , 219)  and  the  file 
minutia  number  2  coordinates  are  (XXF , YYF) = (188 , 211) .     Now,  referring 
to  Figure  21,  suppose  that  the  search  minutiae  are  superimposed  on  the 
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file  minutiae  so  that  search  mlnutla  2  sits  on  top  of  file  minutia  2. 
This  is  illustrated  in  Figure  26.     The  next  ranking  pair  of  minutiae  in 
Figure  25  are  search  4  on  file  5.     Figure  26  shows  the  X  and  Y 
differences  XI  and  YI  as  they  are  computed  in  box  58.006.     The  dis- 
placement of  search  and  file  minutiae  are  shown  as  XIJ  and  YIJ.  Since 
|xi|<|yi|,  control  goes  to  box  58.012  where  the  specific  displacements 
SYT=-.033  and  WXT=.217  are  computed.     If  the  conditions  of  boxes  58.013 
and  58.014  are  satisfied,  these  are  suiraned  as  SY  and  WX  and  a  count  is 
kept  by  ISY.     The  difference  between  the  angles  of  minutiae  search  4 
and  file  5,  TI,   is  summed  as  SA  and  counted  as  LA  in  boxes  58.018  and 
58.019.     (Subroutine  ADA  is  shown  in  Flow  Chart  19.) 

Control  then  Increases  Kl  and  moves  to  box  58.003  where  the  next  pair 
of  minutiae  in  Figure  25,  3  on  3  is  treated.     See  Figure  27.  Since 
XI|>|yi|   control  goes  to  box  58.008  where  SXT=-.186  and  WYT=-.100 
and  these  are  summed  with  similar  values  as  SX  and  WY  in  box  58.011. 
The  last  pass  through  box  58.003  where  (i,  j)=(l,  1)  produces  another 
trip  through  box  58.012.     The  values  of  the  various  variables,  as 
control  makes  three  passes  through  point  AA  in  Flow  Chart  23, are  shown 
in  Figure  28.     When  Kl  reaches  4,  box  58,007  sends  control  up  the  left 
side  of  the  flow  chart  where  the  tensor  components  GOV,  SIW,  COW,  SIV 
and  the  rotation  angle  NUU  are  computed.     For  example: 

C0V=1  -  =  1.186. 

COV  is  the  amount  by  which  one  print  is  stretched  in  the  X  direction 
with  respect  to  the  other  print  per  unit  of  length  in  the  X  direction. 
SIW  is  the  amount  by  which  one  print  is  stretched  in  the  Y  direction 
with  respect  to  the  other  print  per  unit  of  length  in  the  X  direction. 
In  other  words,  COV  and  SIW  deal  with  distortions  in  the  predominately 
horizontal  area  shown  in  Figure  29.     In  a  similar  way  COW  and  SIV  deal 
with  the  predominately  vertical  area  with  respect  to  the  search 
minutia  at   (XXU,  YYU) . 

The  angle  NUU  is  the  mean  of  the  differences  between  the  various  pairs 
of  minutiae  angles  that  were  used  in  computing  the  tensor.  Control 
then  moves  via ^ to  Flow  Chart  24.     The  values  of  the  tensor  being 
illustrated  are  shown  in  Figure  30. 

5.     Transformation  of  Search  Minutiae 

The  search  minutiae  will  now  be  moved  to  new  locations  as  directed  by 
the  tensor  and  some  will  be  eliminated.     The  new  X,  Y,  0  values  will 
then  be  used  to  compute  a  new  S  table  of  DX,  DY,  D9  values  and  control 
will  return  to  Flow  Chart  21  for  computation  of  the  final  matching 
score . 

Box  61  in  Flow  Chart  24  computes  the  XA  and  YA  distance  from  the 
reference  search  minutia  at   (XXU,  YYU)  to  each  of  the  other  search 
minutiae.     If   |xa|>_|ya|   box  61.2  computes  a  new  position  (X,  Y)  for 
this  other  search  minutia,   i,  utilizing  the  tensor  components  COV  and 
SIW.     If   [xa|<|ya|   box  61.3  does  the  job  using  COW  and  SIV.     A  new 
angle  TH  is  computed  for  this  minutia  by  boxes  61.4  -  65.     The  location 


(X,  Y)   is  where  the  search  minutia  i  would  stand  in  the  file  minutiae 
space.     Box  66  then  constructs  a  small  box  about  this   (X,  Y)  location 
as  shown  in  Figure  31.     Boxes  67  -  68  then  search  the  appropriate  area 
of  the  file  minutiae  table  for  any  minutia  j  that  falls  within  this 
box.     If  the  difference  between  the  angles  of  these  prospective  mating 
minutia  is  less  than  ANK  in  box  78,   then  box  79.1  enters  their  DX,  DY, 
D0  values  in  the  new  S  table.     In  practical  applications  this  table 
tends  to  be  limited  to  true  mating  minutiae. 

To  illustrate  the  above  transformation,  consider  the  first  search 
minutiae  in  Figure  22.     Box  61  computes  XA=181-168=13  and 
YA=219-251=-32.     Since   |xa|<|ya|,  box  61.3  computes  the  new  location 
for  this  search  minutia  as:     X=188+(-32) ( . 233)-13=167 . 5  and 
Y=211+32xl . 048=244 . 5 .     This  places  the  search  minutia  almost  perfectly 
on  top  of  the  file  minutia.     The  relocated  search  minutiae  coordinates 
are  shown  in  Figure  32  along  with  the  resulting  S  table.     The  super- 
imposed minutiae  of  Figure  21  after  the  tensor  transofmation  now  appear 
as  shown  in  Figure  33.     After  resetting  the  parameter  KR  to  a  much 
smaller  BKR  =  5,  control  returns  via  B  to  Flow  Chart  21. 

—  6.     Computation  of  Matcher  Score  RS 

The  computer  proceeds  through  Flow  Chart  21  obtaining  the  following 
values  for  the  DX  and  DY  table  shown  in  Figure  32. 


i 

,1 

TR 

KR-TR 

R 

2 

1 

1 

4 

4 

3 

1 

1 

4 

8 
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0 

5 

13 

4 

1 

1 

4 

17 

2 

2 

3 

20 

3 

2 

3 

23 

Control  then  passes  via  C  to  Flow  Chart  22,  and,  since  KSW  now  equals 

one,  goes  to  box  34.03  where  R  is  rounded  to  IIS  for  possible  future 

use.     The  final  RS  score  is  computed  in  box  34.04.     S0LD=4  was  the 

value  of  S  from  the  previous  pass  through  Flow  Chart  21  and  was  saved 

in  box  34.02.     IMXB=4  and  IMXC=5  are  the  numbers  of  minutiae  in  the 

search  and  file  sets.     S=4  is  the  present  number  of  entries  in  the  S 

uT  ^  Rq  =  8500x23x4  _  qiis 

table.     Therefore:  .    ,    ^        -  y//3. 

4x4x5 

7 .  Conclusion 

In  searching  50  cards  against  a  file  of  100  cards  all  of  which  are  from 
the  ulnar  loop  section  of  file,   the  M82  made  a  significant  improvement 
in  the  hit  rate  when  compared  with  the  M40  matcher.     The  M40  software 
matcher  in  the  UNIVAC  1108  computer  matches  a  pair  of  prints-  in  about 
2  milliseconds  compared  with  about  150  milliseconds  required  by  the 
M82  matcher.     The  M82,  like  its  M40  predecessor,  can  be  built  in  high- 
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speed  integrated  circuit  hardware  with  a  significant  improvement  in 
speed  of  performance. 

The  computer  programs  used  in  this  work  were  written  and  managed  by  J. 
F.  Rafferty  of  the  National  Bureau  of  Standards.     The  author  is  indebted 
to  R.  M.   Stock  of  the  Federal  Bureau  of  Investigation  for  advice  and 
assistance.     The  work  was  financially  supported  by  the  FBI. 
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Figure  8.     Portion  of  an  Arch  Pattern  Near  Center  Point  "+" 
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Figure  13. 
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Figure  21.     Search  Minutiae  •  Superimposed  on  File  Finutiae  ^ 
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Figure  24.     The  S  Table 
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Figure  28.     Tensor  Development 
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